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►j,. ■ We estimate in the QCD sum rule approach the amount of the nucleon 

spin carried by the gluon angluar momentum: the sum of the gluon spin and 
L-. , orbital angular momenta. The result indicates that gluons contribute at least 

CN ' one half of the nucleon spin at scale of 1 GeV 2 . 
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Ever since the publication of the EMC measurement on the fraction of the nucleon spin 
carried by the quark spin |I[], there has been a tremendous activity in the field of the spin 
Qh! structure of the nucleon . One of the central questions is how the spin of the nucleon is 

distributed among its constituents || . After much debate, many agree now that a substential 
fraction of the nucleon spin comes from sources other than the quark spin, i.e., quark orbital 
and gluon angular momenta. Recently, several proposals have been made in the literature 
to measure the amount of the spin carried by the gluon helicity AG |4|] . 

In this Letter, we present a QCD sum rule calculation || of the amount of the nucleon 
spin carried by gluons, or equivalently by quarks because, by definition, their sum is 1/2. 
Our calculation is motivated by the possibility of measuring these quantities through deeply- 
virtual Compton scattering proposed by one of us |J. The method we use has been applied 
successfully to calculate a similar quantity — fractions of the nucleon momentum carried by 
quarks and gluons 0,§. Our result shows that the gluon angular momentum, the sum 
of gluon helicity and orbitial angular momentum, contributes at least 50% of the nucleon 
spin, suggesting that the nucleon contains nontrivial gluon configurations carrying nonzero 
angular momentum. 

The angular momentum operator in QCD can be written in an explicitly gauge invariant 
form 0, 

</qcd = J d 3 x [ -^775^ + ^\x x (-iD))ip 

+ xx (Ex B)] . (1) 
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where flavor and color indices are implicit. The first term can be interpreted as the quark 
spin contribution, although its matrix element is actually the singlet axial charge. The 
second term, where the covariant derivative is D = d + igA, is the canonical orbital angular 
momentum of quarks. The word "canonical" stems from the canonical momentum for quarks 
in a background gauge field. The last term is the total angular momentum of gluons, as is 
clear from the appearance of the Poynting vector. [In pure gauge theory without quarks, this 
term generates the spin quantum numbers for glueballs]. According to the above expression, 
we can write down a gauge-invariant spin sum rule for the nucleon, 
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A£(jj?) + L q tf) + J 9 tf) , (2) 



where /z 2 is a scale at which the operators are renormalized, or more physically the nucleon 
wave function is probed. The first term is what has been measured in polarized deep- 
inelastic scattering [!],§]. The second and third terms represent quark orbital and gluon 
contributions, respectively. We also introduce the notion of the total quark contribution, 
J q = AS/2 + L q , the sum of spin and orbital. By definition, both J q (fi 2 ) and J 9 (/i 2 ) are 
gauge-invariant if gauge-invariant regularization and renormalization schemes are used. In 
the light-like gauge A + = 0, J 9 (/i 2 ) can be written as a sum of the gluon helicity AG(// 2 ), 
measurable in polarized high-energy scattering 0], the gluon orbital angular momentum, as 
well as a term from quark-gluon interactions ||. 

Before formulating the sum rule calculation, it is instructive to review a derivation of 
Eq. (1). The angular momentum operators of QCD are identified with the generators of the 
Lorentz group: J'"', which in turn are defined from the angular momentum density M^ va 
through, 

Jf» = fd 3 xM°^{x) . (3) 



The angular momentum density can be expressed in terms of the symmetric, conserved 
energy-momentum tensor T a @, 

M^ a = T» a x u - T^x a . (4) 

The energy-momentum tensor of QCD can be written as a sum of the quark and gluon parts, 

T af3 = T a0 + T a/3 = Ifryfy Q0) ^ + flg^F 2 - F^F^) , (5) 

where (a/3) means symmetrization of the indices. It is then simple to see that the quark and 
gluon parts of the angular momentum operators in Eq. (1) are derived from Eqs. (H) and (|J) 
by substituting in the quark and gluon parts of the energy-momentum tensor, respectively. 
According to the above, we can formulate the sum rule calculation of Jg^fi 2 ), or equiva- 
lently J q (fi 2 ), in terms of the energy- momentum tensor T^g. Consider the following three- 
point correlation function in the QCD vacuum, 

W% va {p) = fd A xd*z (0\T[n{x)f]{0)M^ a {z)}\0) e ip ' x , (6) 

where M^ va is defined as in Eq. (|j) with T a/3 replaced by its gluonic part, and n(x) is the 
interpolating field for the nucleon, which we choose to be ||11||, 



v (x) = e ijk (u iT C-f a ui) 757a rf fc . (7) 

W pua contains a nucleon double-pole contribution, with its residue proportional to J 5 (/^ 2 ), 



(p 2 - m 2 N f 



w r = ( j_y^ 2 ^ u ^ a )+-, (8) 



where ellipses include nucleon double poles of different Dirac structures, nucleon single poles, 
and other dispersive contributions. A^ is the nucleon decay constant corresponding to the 
interpolating current, 

(0\ V (0)\N(p)) = X N U(p) . (9) 

In the following, we first calculate W pua in the deep-Euclidean region — p 2 » Aq CD us- 
ing operator product expansion (OPE), from which we attempt to extract the double-pole 
residue J 9 (/i 2 ). 

To ensure J g ([i 2 ) + J g (^ 2 ) = 1/2 in the sum rule calculation, we use an implicit form of 
Ward identity, 

T af3 = d p {T pp x a ) - x a d p T p/3 , (10) 

to we rewrite the three-point correlation function as, 

W^ a = f d 4 xd 4 z (0\Tr](x)f](0)z p [z a d p T pv - z v d p T pa ) |0)e ip;c . (11) 

From Eq. ([|), we find, 

d p T pv = -$gF va y a i/> + ... , (12) 

where ellipses denote terms vanishing after using gluon's equations of motion. Thus, we 
arrive at a new form of the Green's function 

W pua = fd^xd^z z^z v (0\Tr]{x)fj{0)d a (z)\0)e ipx - {v <-► a) + ... , (13) 

where O a (z) = ipgF al3 'jpilj(z). If one goes through a similar procedure for a correlator with 
the quark part of the energy momentum tensor, one finds that it can be reduced to the same 
term with a negetive sign plus a two-point nucleon correlation function with a double-pole 
residue 1/2. 





(a) (b) 

FIG. 1. Perturbative diagrams. Dashed line denotes gluon. (Permutations are not shown) 



The Green's function in the deep Euclidean space can be calculated in OPE because 
of asymptotic freedom. The first term in such an expansion is the usual perturbative con- 
tribution, which is infrared finite due to the finite external momentum p 2 . There are two 
perturbative diagrams as shown in Fig. 1. We find the contribution from the first diagram 
as, 

-J—^^--\n^] P \ (14) 

vr 5 ^144 /i 2 36 p 2 j y ' V ; 

where and henceforth we omit the structure factor 2ip p/ j l/ p 1 r y a . A calculation for the second 
diagram ("salboat") is rather tedious. Since in the final result the (typical) contribution 
from the first diagram is small (less than 10%), we discard this "sailboat" contribution in 
the following study. 






(a) (b) (c) 

FIG. 2. Dimension-4 power corrections: local (a,b) and bilocal (c). Shaded circles mark vacuum 
fields 

The next term in OPE comes from dimension-four vacuum condensates. Diagrams from 
Fig. |2|a,b are found to contribute, 

1 : (V>(4 + I + ta ^, (15) 



1447T 2 p 2 7T y p? 6 — q 2 / 

where q 2 is an infrared regulator which represents the momentum flow through the operator 
O a . The infrared logarithm arises from large separations of point z from and x. To take 
into account the contribution in this region properly, one must first expand the product of 
the interpolating current, 

T V (x)fj(0)='£C n (x)d n , (16) 



(where O n are a set of local operators) resulting in so-called bilocal power corrections (12 
(see Fig. ^p). The relavant local operator in this case is a dimension-five one, 

6l pp ' = 2ugF x[f, Y' ] u - 2id [p (u L? l p ' ] u) + u $5* a^'u 

+ua pp ' D X l) u + -ugFPP'^u + -dgF pp '^ x d , (17) 

where [pp'] denotes antisymmetrization of the two indices. The operator yields a contribution 
to W^ a , 

n (0,/i 2 ), (18) 



12tt 2 p 2 
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where n (g 2 ,/i 2 ) is a two-point correlation function between O b pp and O a , and p 2 is an 
ultraviolet regulator to be defined below. 

To calculate n (0, /x 2 ), we again use the sum rule approach. We first work out an 
operator-product expansion for n (g 2 ) in the deep Euclidean space, 



n (gV 2 ) 



q In — - + — (— F )ln 
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q* yg" I92ir 2 q 2 

On the other hand, we write a dispersion integral for Ho^ 2 ,/! 2 ) valid for all q 2 |14 

1 f^ (JS p(s) - p P ert{s) 

f Jo s 



(g"G") + ... (19) 



n (g 2 ,/i 2 



s — q 
where the upper limit defines the ultraviolet cut-off and 

/ X Ot s q 

Ppert(, s J 



607T 2 



(20) 



(21) 



Ilo(g 2 ) defined in this way vanishes in perturbation theory and its first power correction 



contributes in the same way as the last term in Eq. (iSf). To find ll (0,/i 2 ), we assume a 
spectral function, 



p{s) = irf R m e R 5(s - m 2 R ) + 9{s - s ) 



a. 



60vr : 



r , 3 + ^F 2 , 
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(22) 



where m R is the mass scale for the exotic 1 h resonace, suspected to lie between 1.3 to 
1.9 GeV |U|. In our estimate, we take m R to be 1.5 MeV. The standard sum rule method 
allows us to extract f R = 1.8 x 10 -3 , which in turn yields Uo(Q,m R ) = 5.0 x lO~ 4 m R . The 
uncertainty of this number is at least a factor of 2 due to unknown m R and the continuum 
threshold sq, which we take to be 1.9 2 GeV 2 . 

The next term in the OPE for W^ ua involves dimension-six vacuum condensates, for 
which we use the factorization assumption. A calculation of the diagrams in Fig. |3|a,b,c 
(and similar ones which are not drawn) give a contribution to W^ va , 
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FIG. 3. Typical local (a,b,c) and bilocal (d) power corrections of dimension 6. 
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where we have kept only logarithmic terms. Small contribution of the first term to the final 
result justifies the approximation. The infrared logarithm in the second term signals that 
the contribution must be replaced by, 



4 Wn ,n ..2 



3p 4 



^(O,^ 2 ), (24) 



where Il 1 (g 2 ,/i 2 ) is a bilocal correlator (see Fig ||d) give involving O a and the dimension- 
seven operator, 

qW = e W e i'?k( D \ u y Cj p u i u j 'Y'Cu i/T + h.c. (25) 

The OPE for n^g 2 ,^ 2 ) at large Euclidean q 2 is, 

„ , Q . 31 a s . . , u 2 mkiuu) , . 

where m 2 , = —(ugF ■ au)/(uu). The higher-order terms in ellipese involve condensates of 
dimension-seven and higher for which we know very little. To get an estimate, we assume 
vector- meson dominace |15[ |, 



IW) = -4^~ 2 • (27) 

m R ~ Q 

Expand the above in q 2 and matching its 1/g 2 term with the OPE in Eq. (|2HD , we find, 

0,(0) = ?f$. . (28) 

We ignore dimension-eight or higher contributions. In the factorization approximation, 
the contributions from dimension-eight condensates (both local and bilocal) are exactly zero. 

Based on the OPE we have developed for W^ va , we attempt an estimate for the J g (m 2 N ). 
The sum rule equation reads like this, 

+ - = r^(T77 ln -^ - ^ In -ir)p - t-— ^(— F ) In -- + 



(m^-p 2 ) 2 '" vr 5V 144 ' /i 2 36 ,u 2 ^ 144ttV n n '\ // 2 6/ 

1.1 xlQ- 3 20a s (^) 2 -p 2 4m 2 (mx) 2 
12ttV 81tt p 4 n /i 2 Om 2 ^ 4 ^ ^ 

Substituting in the standard values for the condenstates at the normalization pont // = lGeV 
(ccf ref. |16| for example): ((a s /n)F 2 ) = 0.012 GeV 4 , (uu) = -0.017 GeV 3 , m 2 = 0.65, 
a,(l GeV) = 0.37, 32vr 4 A 2 v = 2.5 GeV 6 , s = 2.25 GeV 2 , we find that the dimension-six 
bilocal term is the dominant contribution. If we keep just this term, we find 

A more careful analysis including other contributions yields, 

J g (l GeV 2 ) = 0.35 ±0.13 . (31) 

where the error reflects the uncertainty of the mass scale in the 1 h channel as well as the 
uncertainty from the sum rule analysis. However, we have no way to know the accuracy of 
the vector meson approximation in estimating the dimension-six bilocal contribution. 



The number we find, J g (l GeV 2 ) ~ 0.35 ± 0.13 or J q {\ GeV 2 ) ~ 0.15 ± 0.13, if taking 
seriously, has interesting implication on the spin structure of the nucleon. It says that 
gluons are at least as important in determining the nucleon spin as quarks, if not more. 
Furthermore, from a recent globle analysis of data on polarized deep-inelastic scattering 
||, one finds the gluon helicity AG(1 GeV 2 ) defined in the infinite momentum frame and 
light-like gauge has a size of 1 to 2 units of angular momentum. If correct, the gluon 
orbital contribution defined in a similar framework must be large and negative and cancel 
a substential part of AG. Such a large cancellation may be caused by the gauge-dependent 
separation of J g into helicity and orbital contributions. On the other hand, one half of the 
singlet-axial charge, or the quark spin contribution, is found to be 0.05lo;o5 0. This leaves 
about 20% of the nucleon spin carried by quark orbital angular momentum. Here no large 
cancellation is present between the quark spin and orbital contributions. 
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